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We examine the claim by Oppenheimer et al. (2001) that the local halo density of 
white dwarfs is an order of magnitude higher than previously thought. As it stands, 
the observational data support the presence of a kincmatically distinct population 
of halo white dwarfs at the >99%confidence level. A maximum-likelihood analysis 
gives a radial velocity dispersionrl of o-^=150^!q kms _1 and an asymmetric drift 

of t>a=176^gQ 2 kms -1 , for a Schwarzschild velocity distribution function with 
<T[7 :0 "v:o"M'=l:2/3:l/2. Halo white dwarfs have a local number density of I.I^q'^ x 
10 — 4 pc — 3 , which amounts to 0.8^q | per cent of the nominal local dark-matter 
halo density and is 5.0^3 2 times (90% C.L.) higher and thus only marginally in 
agreement with previous estimates. We discuss several direct consequences of this 
white-dwarf population (e.g. microlensing) and postulate a potential mechanism to 
eject young white dwarfs from the disc to the halo, through the orbital instabilities 
in triple or multiple stellar systems. 



1 Introduction 

Recently, Oppenheimer et al. (2001; O01 hereafter) found 38 white dwarfs 
(WD) in a sample of 99 high proper motion WDs, which were claimed to have 
kinematics inconsistent with that of both the old stellar disc and the thick disc, 
and therefore form a newly discovered halo population with an inferred density 
at least an order of magnitude higher than previously thought (e.g. Gould et 
al. 1998). Reid et al. (2001) have criticized this result and conclude that 
these WDs form the high-velocity tail of the thick disc, based on a comparison 
of their U-V velocity distribution with that of M-dwarf stars. Unfortunately, 
neither conclusion is statistically supported. 

In this proceeding, we discuss a maximum-likelihood analysis of the com- 
plete sample of 99 WDs and derive the local phase-space density of thick-disc 
and halo WDs. This then allows us to draw more robust conclusions. For a 
discussion of the selection of the sample, we refer to O01. Details of the likeli- 
hood analysis and phase-space density estimate can be found in Koopmans & 
Blandford (2001; KB01 hereafter). For clarity, throughout the proceeding the 
low- and high-velocity dispersion components are referred to as the thick-disc 
and halo population, respectively. This does not imply that all WDs could not 
have originated from the thin or thick disc (e.g. Hansen 2001; KB01). 

"Hereafter, if not otherwise indicated, errors give the 90% statistical confidence level. 
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2 The White Dwarf Velocity Distribution Function 

We model the local (x*o) velocity distribution function (VDF) as a superposition 
of two Schwarzschild (i.e. Maxwellian) VDFs, i.e. /b(xo,u). We assume a con- 
stant space density of WDs throughout the surveyed volume. Each component 
has two free parameters, the radial velocity dispersion (au) and the asymmet- 
ric drift (v a ). In addition, the parameter r„ is the ratio of the thick-disc to halo 
WD number densities (n). We assume that the vertex deviations are zero and 
that the ratios of the radial ([/), azimuthal (V) and vertical (W) velocity dis- 
persions of the ellipsoidal VDFs are c7[/:ov:<7vk=1:2/3:1/2, in agreement with 
observations (see KB01 for references). The probability that an observed WD 
with a velocity vector on the sky, p = (vi,Vb), is drawn from this VDF becomes 
P = C( v l + v 2 b ) 3/2 .fB(v = p + v r f)dv r / f v (v? + vlf/ 2 f B {v = p + v r r)d*v, 
where v r is the radial velocity along the line-of-sight and vi and Vb are the veloc- 
ities projected on the sky in Galactic coordinates. This assumes that the WDs 
are all proper-motion limited (^90% arc). V is then the velocity space in which 
the WDs could have been found, given the restriction that vf + v\ > (/j, r) 2 , 
where r is the distance to the WD and juo is the survey's lower limit on the 
proper motion (/io ^0.33 "yr _1 ). For the ~10% magnitude limited WDs, we 
find that their maximum detection volumes (Vmax) are nearly identical. For 
these, we modify the likelihood function, although the differences in param- 
eter estimates are within the errors quoted below. By varying the five free 
parameters and optimizing the log-likelihood, log(Pj), of the sample, 

we solve for their most likely values and their error range. The results (exclud- 
ing 3 thin-disc WDs) are: (a) CT[/ i d=62+ 10 kms -1 , (b) 0-[/ ; h=15O+ 4O kms -1 , 
(c) v a4 =50tu kins" 1 , (d) v a ,h=176±^ 2 km s" 1 and (e) r n =16±?° (KB01). 

3 The Local Halo White Dwarf Density 

To normalize the local phase-space density of halo plus thick-disc WDs, i.e. 
^owd ' /b(xo,w), it remains to estimate their local density, Iqwd - 

We do 

this in the conventional way, by summing the 1/Knax values of the halo plus 
thick-disc WDs, where V max is the smallest of the two volumes in which a WD 
could have been detected, when limited either by its proper motion or by its 
magnitude (e.g. O01; KB01). We find that ~90% of the WDs arc proper- 
motion limited and three WDs are very likely associated with the thin disc. 
The latter WDs are removed from the sample, as we are only interested in 
the halo plus thick-disc density, for which we determined the VDF. We find 
n* d +£ = (1.9 ± 0.5) x 1CT 3 pc- 3 , of which n^ WD = l.lto.7 x 10_4 P c ~ 3 
belongs to the halo, given r n found from the maximum-likelihood analysis of 
the sample. The local halo plus thick-disc density estimate agrees within ~1 a 
with the estimate from Reid et al. (2001). The halo density is 5.0l|;2 times 
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higher than the best previous estimate of h§ WD = 2.2 x 1CP 5 pc~ 3 from Gould 
et al. (1998), inferred from the mass function of halo stars. Locally, our density 
estimate amounts to 0.8^q'5% of the nominal local dark-matter halo density 
of 8 x lfr 3 Mq pc" 3 (Gates et al. 1995). 

4 The Global Halo White-Dwarf Density 

Given the local phase-space density of halo plus thick-disc WDs and a potential 
model for the Galaxy, we can estimate the global white-dwarf density in a large 
part of the Galactic halo, using Jean's theorem (e.g. May & Binney 1986). For 
simplicity, we assume a spherical logarithmic potential, <I>(r) = v^\n(r/r c ), 
with rotation velocity v c =220 kms -1 and a solar radius of r c —8 kpc. For the 
halo WDs, which are less affected by the disc/bulge potentials this is probably 
a reasonable assumption, although this does not hold for the lower-velocity 
thick-disc WDs. A analytic expression for the phase-space density of halo 
WDs can then be derived for any point in halo that is in dynamic contact with 
the Solar neighborhood (KB01) . From this phase-space density model, we find 
an oblate (q = (c/a) p ~ 0.9) distribution of halo WDs with a total mass inside 
50 kpc of ~2.6xl0 9 M Q and a radial mass profile n(r) oc r~ 3 - . We expect q 
to decrease further if a proper flattened Galactic potential is used. This halo 
WD mass amounts to ~0.4% of the total Galactic mass inside 50 kpc, ~4% of 
the Galactic stellar mass, or f2 WD <~ 10~ 4 . Including disc WDs, the total mass 
of WDs in the Galaxy is ^9% of the stellar mass, in agreement with standard 
population synthesis models (Hansen 2001; KB01). 

5 The LMC Microlensing Optical Depth 

Given our analytic expression of the halo WD phase-space density, we can 
integrate over the line-of-sight towards the Large Magellanic Cloud (LMC) to 
estimate the microlensing optical depth of this population. We find r^ D ~ 
1.3 x 10~ 9 , which is ~ 10 2 times lower than observed (Alcock et al. 2000). 
The thick-disc WDs have a ~3 times higher optical depth. The integrated halo 
plus thick-disc WD optical depth is therefore still 1-2 orders of magnitude too 
small. 

6 Conclusions, Discussion & Future 

Modeling of the local phase-space density of halo plus thick-disc WDs indicates 
that there is a population of pressure-supported halo WDs with a local density 
5.0^32 times higher than previously estimated, although globally (r < 50 kpc) 
the mass contribution is negligible (<~0.4%). Both the low- and high- velocity 
populations have similar color, magnitude and ages distributions (Hansen 2001; 
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KB01). Two particularly interesting surveys to search for more high proper 
motion WDs are: (i) a wide survey with the Advanced Camera for Surveys 
(ACS) on the Hubble Space Telescope, which we estimate should find ~5 WDs 
(KB01) and (ii) a similar survey as that by O01 towards the Galactic anti- 
center to constrain the WD velocities perpendicular to the Galactic plane. 

The white-dwarf data set from O01 and results derived from it (i.e. our 
most likely model) have met a number of interal and external consistency 
checks (KB01) and agree with all observational constraints (e.g. direct obser- 
vations, metal pollution of the ISM, microlensing, etc.) known to us. The 
color-color relation of the WDs shows them to be quite young (Hansen 2001) 
and have a birth rate roughly that expected from Galactic population synthesis 
models (KB01). Why then do a large fraction of these WDs have such high 
spatial velocities? To explain this we postulate the following; the majority of 
high-velocity WDs in the Galactic halo originate from multiple (i.e. N > 3) 
stellar systems in the Galactic disc. In case of a striple stellar system, stellar 
evolution (i.e. mass loss and/or transfer) of the (probably more massive) inner 
binary stars changes the ratios of orbital periods and eventually destabilizes 
the system. At that point, the lightest star, a WD in the inner binary, is 
ejected from the system into the halo - possibly through a slingshot of the 
outer star - leaving behind a recoiled compact binary (KB01 ). This postulate 
implies a direct relation between the radial surface density of the stellar disc, 
the star-formation rate and population synthesis models, which all seem to be 
in agreement (KB01), but still need considerable study. 
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